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ABSTRACT 


Reaction turbines basically Francis turbines and propeller/Kaplan turbines are suitable for medium and 
low head hydropower sites. The management of the small hydropower plants is an important factor, for 
achieving higher efficiency of hydro turbines with time. Turbines show declined performance after few 
years of operation, as they get severely damaged due to various reasons. One of the important reasons is 
erosive wear of the turbines due to cavitation. Reaction turbines, however are more prone to cavitation 
especially Francis turbines where a zone in the operating range is seriously affected by cavitation and 
considered as forbidden zone. Cavitation is a phenomenon which manifests itself in the pitting of the 
metallic surfaces of turbine parts because of the formation of cavities. In the present paper, studies 
undertaken in this field by several investigators have been discussed extensively. Based on literature 
survey various aspects related to cavitation in hydro turbines, different causes for the declined 
performance and efficiency of the hydro turbines and suitable remedial measures suggested by various 


investigators have been discussed. 


© 2009 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Energy is the basic requirement for economic development. 
Every sector of a country’s economy (industry, agricultural, 
transport, commercial and domestic) needs input of energy. The 
installed electricity generating capacity in India at presentis nearly 
128 GW which is generated 66% through thermal, 26% through 
hydro, 3% through nuclear and 5% through other renewables. 
According to International Energy Agency (IEA), a threefold rise in 
India’s generation capacity is expected by 2020 [1]. As the non- 
renewable fossil energy sources continues to deplete, and realizing 
the summits held at Brazil and Kyoto, to reduce the greenhouse 
gas emissions, hydropower has moved towards the top power 


* Corresponding author. Tel.: +91 1332 285841; fax: +91 1332 273517/60. 
E-mail address: rajsafah@iitr.ernet.in (R.P. Saini). 


1364-0321/$ - see front matter © 2009 Elsevier Ltd. All rights reserved. 
doi:10.1016/j.rser.2009.07.024 


development option to meet the increasing energy demand. Small 
hydropower (SHP) has grown substantially in the last 10 years. 
Among all the renewable energy sources available, small hydro- 
power is considered as the most promising source of energy. The 
hydraulic turbine being the heart of a SHP plant transforms the 
potential energy of water into mechanical energy in the form of 
rotation of shaft. Hydraulic turbines can be broadly classified into 
two categories according to action of water on moving blades. 
Hydro turbines classification based on action of water over runner 
are impulse turbines and reaction turbines. Impulse turbines are 
high head turbines and operate at atmospheric pressure. At some 
of hydropower sites where head is relatively low and discharge is 
of medium or high order it becomes difficult to operate an impulse 
type turbine as sufficient speed may not be obtained to drive the 
end use machine at such sites. It is therefore reaction turbines are 
suggested to be used and these considered high speed turbines 
under medium and low head. These turbines are closed turbines 
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operate under variable pressure. Further these turbines are more 
prone to cavitation and a zone in the operating range is seriously 
affected by cavitation and considered as forbidden zone especially 
in case of Francis turbine. The management of the small hydro- 
power plants for achieving higher efficiency of hydro turbines with 
time is an important factor, but the turbines show declined 
performance after few years of operation as they get severely 
damaged due to various reasons. One of the important reasons is 
erosive wear of the turbines due to high content of abrasive 
material during monsoon and cavitation [2]. According to the 
Bernoulli’s equation, if the velocity of flow increases, the pressure 
will fall. In case of liquid, the pressure cannot fall below vapor 
pressure, the vapor pressure is defined as the pressure at which 
liquid is vaporized at a given temperature, which depends upon the 
temperature and height above mean sea level of the site. Whenever 
the pressure in any turbine part drops below the evaporation 
pressure, the liquid boils and large number of small bubbles of 
vapors are formed. It may happen that a stream of water cuts short 
of its path giving rise to eddies and vortices which may contain 
voids or bubbles. These bubbles mainly formed on account of low 
pressure, are carried by the stream to higher pressure zones where 
the vapors condense and the bubbles suddenly collapse, as the 
vapors are condensed to liquid again. This results in the formation 
of a cavity and the surrounding liquid rushes to fill it. The streams 
of liquid coming from all directions collide at the center of cavity 
giving rise to a very high local pressure whose magnitude may be 
as high as 7000 atm. Formation of cavity and high pressure are 
repeated many thousand times a second. This causes pitting on the 
metallic surface of runner blades or draft tube. The material then 
fails by fatigue, added perhaps by corrosion. Some parts of turbine 
such as runner blades may be torn away completely by this 
process. The phenomenon which manifests itself in the pitting of 
the metallic surfaces of turbine parts is known as cavitation 
because of the formation of cavities. Due to sudden collapsing of 
the bubbles on the metallic surface, high pressure is produced and 
metallic surfaces are subjected to high local stresses. Thus the 
surfaces are damaged [3]. Prof. D. Thoma suggested a dimension- 
less number, called as Thoma’s cavitation factor o (sigma), which 
can be used for determining the region where cavitation take place 
in reaction turbines: 


_ Hyp- Hs 
=e (1) 


where H, is the barometric pressure head in m of water, H, is the 
suction pressure at the outlet of reaction turbine in m of water or 
height of the turbine runner above the tail water surface, H is the 
net head on the turbine in m of water [4]. 

The value of Thoma’s cavitation factor (o) for a particular type of 
turbine is calculated from the Eq. (1). This value of Thoma’s 
cavitation factor (o) is compared with critical cavitation factor (o<) 
for that type of turbine. If the value of ø is greater than o, cavitation 
will not occur in that turbine. 

For safe operation (cavitation-free) of turbine, it is evident that: 


o>Te (2) 


o 


The following empirical relationships are used for obtaining the 
value of T. for different turbines [5]. 

For Francis turbines, 
Te = 431 x 1078 x N2 (3) 


For propeller turbines, 


1 


Te = 0.28 4 Z 
7.5 x (Ns/380.78) 


where N, is the specific speed of turbine. 
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Fig. 1. Variation of efficiency with respect to cavitation factor, o. 


Fig. 1 shows the variation of Thoma’s cavitation factor (o) with 
efficiency. With increase in suction head (H,), a corresponding 
decrease in o is obtained and this has no effect on the turbine 
performance as is seen from constant efficiency trend. A stage is 
however reached when any further increase in suction head (Hs) 
deteriorates the turbine performance and the turbine efficiency 
falls. The critical value of cavitation coefficient, Te is determined by 
this point marking the change in the trend of efficiency [5]. 

The present paper discusses the studies carried out by various 
investigators in order to determine the effect of cavitation and 
identify the gaps for future studies. 


2. Theoretical investigations 


Kjolle [6] studied the causes of damages in hydro turbines and 
found that the main causes of damage of water turbines were due 
to cavitation problems, sand erosion, material defects and fatigue. 
The turbine parts exposed to cavitation are the runners and draft 
tube cones for the Francis and Kaplan turbines. The effect of 
cavitation erosion was found to be reduced by improving the 
hydraulic design and production of components, adopting erosion 
resistant materials and arrangement of the turbines for operations 
within the good range of acceptable cavitation conditions. 

Duncan [7] provided information about cavitation, cavitation 
repair, cavitation damage inspection, date of inspection, cause of 
pitting, runner modifications, cavitation pitting locations, methods 
of cavitation pitting repair, and areas of high stress in runner. 
Typical areas of cavitation pitting were found as shown in Fig. 2. 

Farhat et al. [8] illustrated the benefits of the cavitation 
monitoring in hydraulic turbines using a vibratory approach. This 
technique was used in a large Francis turbine in order to validate as 
light modification of its distributor, which was intended to reduce 
the cavitation aggressiveness and thereby the related erosion. 
Cavitation-induced vibrations were measured in fixed parts of the 
turbine prototype and compared to those measured in a similar 
and non-rehabilitated turbine. 

Karimi and Avellan [9] presented a new cavitation erosion 
device producing vortex cavitation. A comparative study between 
various cavitation erosion situations was carried out to verify the 
ability of this vortex cavitation generator to produce realistic 
cavitation erosion with respect to that observed in hydraulic 
machinery. Hardened superficial layers in specimens exposed to 
flow cavitation were found a thicker than those in vibratory 
cavitation, which leads to higher erosion rates. 

Huixuan et al. [10] carried out an on-line monitoring system 
for water turbines. Both audible sound (20Hz-20kHz) and 
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Fig. 2. Francis turbine—typical areas of cavitation pitting. 


ultrasound (50-300 kHz) were continuously monitored. The signal 
characteristics such as standard deviation, noise level and 
frequency compositions were evaluated. The evaluation results 
were stored in database in association with the operating condition 
donated by water head and wicket gate opening (or power output). 
According to the frequency characteristics, the sound emitted by 
cavitation was distinguished from the others such as water flow 
sound and mechanical sound. The cavitation intensity on different 
water head and power were traced out. The degree of cavitation 
erosion was estimated according to the cavitation intensity at a 
fixed operating condition such as rated power and designed water 
head. 

Alligné et al. [11] studied the hydroelectric power plants for 
their ability to cover variations of the consumption in electrical 
power networks. In order to follow this changing demand, 
hydraulic machines were subjected to off-design operation. In 
that case, the swirling flow leaving the runner of a Francis turbine 
may act under given conditions as an excitation source for the 
whole hydraulic system. In high load operating conditions, vortex 
rope behaved as an internal energy source which leads to the self- 
excitation of the system. The aim of this paper was to identify the 
influence of the full load excitation source location with respect to 
the eigen modes shapes on the system stability. 


3. Experimental studies 


Escalera et al. [12] carried out an experimental investigation in 
order to evaluate the detection of cavitation in actual hydraulic 
turbines. The methodology was based on the analysis of structural 
vibrations, acoustic emissions and hydrodynamic pressures 
measured in the machine. The results obtained for the various 
types of cavitation found in the selected machines were presented 
and discussed in detail. Various types of cavitation in Francis 
turbines were found are as shown in Fig. 3. In traveling bubble, the 
generalized Rayleigh—Plesset equation was found valid approx- 
imation of the bubble growth and it can be solved to find the radius 
of the bubble, R(t); provided that the bubble pressure, P(t); and 
the infinite domain pressure, P,.(t); are known. The equation is 
expressed as: 


dRg 3 (/dRg\* 4vdRg 2 
B ( 2) i B, Y (5) 


Pa(t) — Palt) _ p ORs | 
p B AE J dt " Rg dt ' (ORB 


where v is the kinematic viscosity, ‘y is the surface tension and 
‘ep’ is the density. 


Bajic [13] carried out noise sampling, signal processing and 
analysis and data processing, analysis and interpretation in vibro- 
acoustic diagnostics of turbine cavitation. These were investigated 
in a series of prototype and model experiments by several weak 
points of the practice and were identified as shown in Fig. 4. 
Improvements and new techniques were developed. These 
techniques enabled extraction of data on cavitation details and 
early detection of detrimental effects met in turbine exploitation. A 
brief review of weak points of the practice, developed improve- 
ments, and new techniques, as well as examples of application was 
presented. Noise power, /(P), at the turbine power P, used as an 
estimate of the cavitation intensity, could be simply modeled by 
the exponential rule: 


M 
(P) = So Im (6) 
m=1 


where Im is its relative amplitude in the fully developed stage. 

Zuo and co-workers [14] carried out a study related to 
cavitation associated pressure fluctuations in hydraulic systems 
which was an important phenomenon that affected the design, 
operation and safety of systems such as hydropower plant, 
pumping station, rocket engine fuel system, water distribution 
network, etc. (Fig. 5). The physical model of cavitation cloud in the 
system was assumed as an exciter, which had its own character- 
istic frequency (fcg,) rather than simply a lumped capacitance (C) 
and, feay was a function of cavitation number, and operating 
condition (Q4, nı). In turbo-machinery system, frequency was 
expressed as: 


cav 7 Seav(O, nm, Qı) (7) 


Wen-quan etal. [15] carried out investigations on turbulent flow in 
a 3D blade passage of a Francis hydro turbine and flow was 
simulated with the Large Eddy Simulation (LES) to investigate the 
spatial and temporal distributions of the turbulence. The 
computed pressures on the pressure and suction sides agreed 
with the measured data for a working test turbine model. In LES, 
large-scale structures could be obtained from the solution of the 
filtered Navier-Stokes equations, in which the structures smaller 
than the grid size need to be modeled is expressed as: 


a ð ag FP ðf (am, 3N] ðt 
at * Ox; E) = — Be, + Gi, KE 7) Ox; 
— 2 j3KUx (8) 


where P= p — w*r?/2, i; is the large-scale relative velocity, p; is 
the pressure divided by the fluid mass density, w is the constant 
angular velocity around the x3-axis, r is the radial distance to the 
x3-axis, &3, is a circular replacement tensor, v is the kinetic 
viscosity, and t;; = ujuj; — uju; is the SGS stresses. 

Escaler et al. [16] carried out work to improve the cavitation 
erosion prediction methodology in hydro turbines by the use of 
onboard vibration measurements taken on the rotating shaft. It 
was discussed whether the use of the vibration from the shaft as an 
alternative to the bearing ones should be an advantage to infer the 
absolute erosive forces taking place on the runner blades. 

Athavale et al. [17] investigated a new “full cavitation model” 
developed for performance predictions of engineering equipment 
under cavitating conditions. All the test cases with cavitation 
showed plausible results (no negative pressures, and good 
convergence characteristics). Computations on the water jet pump 
for non-condensable gas concentrations showed sizeable changes 
in the pump head developed. The vapor transport equation 
governing the vapor mass fraction, f, is given as: 


2 (pf) + VOV f) = YY f) + Re — Re (9) 
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Fig. 3. Different types of cavitation in Francis turbines: (a) leading edge cavitation, (b) traveling bubble cavitation, (c) draft tube swirl and (d) interblade vortex cavitation. 


where p is the density, V is the velocity vector, lis the effective 
exchange coefficient, and Re and R. are the vapor generation and 
condensation rate terms. 

Nicolet et al. [18] carried out an investigation on scale model 
of high specific speed Francis turbines. The analysis of the 


Fig. 4. The sensors placed on the 20 guide vanes react cavitation in various locations 
around the spiral. 


resulting pressure fluctuation in the entire test rig shows 
significant pressure amplitude mainly at 2.46 fn, which evidences 
the excitation mechanism. A component of pressure fluctuation at 
~2.5 fn frequency was identified all along the draft tube walls, the 
source of those pressure fluctuations being located at the inner 
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Fig. 5. Cavitation resonance from the UM Venturi. 
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Fig. 6. Cumulative weight loss of various industrial alloys. 


part of the draft tube elbow. The analysis of the pressure 
fluctuations phases for the ~2.5 fn frequency provides a way to 
determine the experimental wave speeds along the draft tube, 
which were the key parameters for a numerical simulation of the 
hydro-acoustic behavior of the test rig. The simulation carried out 
for the full test rig, taking into account piping, circulating pumps 
and the scale turbine model with the elbow draft tube showed that 
the ~2.5 fn frequency value corresponds to an eigen frequency of 
the system. 

Hart and Whale [19] carried out work on improved weld 
surfacing alloy which was developed and tested to resist 
cavitation-erosion in hydro turbines. Typical wear characteristics 
experienced in laboratory testing correlated to actual service 
conditions. A metallurgical evaluation showed that a high 
strain, work hardening austenitic stainless steel produces superior 
resistance to cavitational erosion. Several industrial alloys were 
evaluated using the vibratory and high velocity cavitation test, to 
produce a new alloy development in weld surfacing. Field testing 
showed an improvement in cavitation—erosion resistance of up to 
800% relative to 308 stainless steel. The cumulative weight loss of 
various industrial alloys are shown in Fig. 6. 

Jean-Francois et al. [20] presented a comprehensive research 
effort aimed at understanding the influence of flow unsteadiness 
on leading edge cavitation. Pressure measurements were made on 
the suction side of an oscillating hydrofoil in the high speed 
cavitation tunnel. Leading edge roughness effects were investi- 
gated for both cases of fixed and oscillating hydrofoils. The cavity 
that detaches on the rough leading edge was thicker while its 
length was approximately the same as for the smooth configura- 
tion. Pressure fluctuations in the cavity closure region were higher 
for the rough configuration. It was observed that the cavity 
detaching on the smooth leading edge at fixed incidence can be 
completely vanish when the hydrofoil oscillates (Fig. 7). 

The term cavitation was used to describe the phenomenon of 
liquid-to-gas and gas-to-liquid phase changes that occur when the 
local fluid dynamic pressures in areas of accelerated flow drop 


Fig. 7. Stable (left) and unstable (right) leading edge cavities. 
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Fig. 8. Typical eroded areas of a Francis runner. 


below the vapor pressure of the local fluid [21]. The liquid-to-gas 
phase change was akin to the boiling of water, except that it 
occurred at ambient temperatures. The gas-to-liquid phase change 
produced extremely high local pressures as vapor cavities implode 
on themselves. Cavitation commonly occurred in hydroelectric 
turbines, generally appearing around guide vanes, wicket gates, 
the turbine runner, and in the draft tube. Usually, cavitation within 
the fluid stream is not damaging to the turbine. However, when 
implosions occur near solid boundaries within the machine, flow 
surfaces can be damaged and eroded. Detection of the cavitation 
phenomenon is straight forward. Large increases in noise, 
particularly in moderately high frequency ranges (15-100 kHz) 
are characteristic of cavitation. In addition, vibration levels 
generally increase. However, in a machine condition monitoring 
program, the simple ability to detect cavitation is not too 
beneficial. The real need is to learn when the cavitation is 
damaging parts of the machine. 

Avellan [22] found out that design, operation and refurbish- 
ment of hydraulic turbines, pumps or pump-turbine were strongly 
related to cavitation flow phenomena, which may occur in either 
the rotating runner-impeller or the stationary parts of the 
machine. The study presented the cavitation phenomena featured 
by fluid machinery including type of cavity development related to 
the specific speed of machines in both pump and turbine mode, the 
influence of the operating conditions, such as load, head and 
submergence. Therefore, for each type of cavitation illustrated by 
flow visualization made at the EPFL testing facilities, the influence 
of cavitation development on machine efficiency, operation and 
integrity were also discussed in Figs. 8-14. 

Grekula and Bark [23] carried out an experimental study on 
cavitation processes in a Kaplan model turbine and studied with 
the aim to identify mechanisms that promote erosive cavitation. 
The studies were carried out with high-speed filming, video filming 
and visual observations with stroboscopic light. A periodic pattern 
of the cavitating tip vortex was observed. Cavitation at blade root 
observed as shown in Fig. 15. The main feature of the pattern was 
that the cavitating vortex was bent towards the blade surface and 
transformed into cloud formations. It was also found that these 


Fig. 9. Typical erosion at the wall of the blade suction side due to inlet edge 
cavitation, shaded area A of Fig. 8. 
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Fig. 10. Typical erosion at the wall of the blade suction side due to inlet edge 
cavitation, shaded area B of Fig. 8. 


cloud formations appears in bands, with a periodicity which 
corresponds approximately to the spacing of the guide vanes. 
Moreover, the tip vortex cavitation was characterized by fine- 
scaled cloud cavitation at the points where it approached the blade 
surface, which indicated an erosion risk. 

Harano et al. [24] investigated two splitter-blade-fitted runners 
installed in the Francis turbines. The splitter-blade-fitted runner 
was a practical application achieved through joint research, 
utilizing the latest fluid-analysis technology, model testing, 
strength analysis, etc., between Hitachi and The Kansai Electric 
Power Co., Inc. The relationship between inlet diameter and outlet 
diameter of Francis turbine is shown in Fig. 16. Results on vibration 
between turbine metal case and draft tube are shown in Fig. 17. 

Szkodo [25] carried out an experimental investigation of 
cavitation resistance of Fe-Cr-Mn coating. The alloy was used in 
mending of machine elements subjected to cavitation. Chromium 
nickel stainless steel OH18N9 was used as the substrate. The 
investigated sample was exposed to cavitation loading at the 
rotating disk facility. The investigations were performed in the 
initial period of the clad damage. The plastic deformation on the 
investigated surfaces was defined using an image analysis. The 
microstructure, chemical composition and phase identification of 
the modified layer were examined using scanning electron 
microscopy (SEM), light microscopy (LM), energy dispersive X- 
ray spectroscopy (EDX) and X-ray diffractometry (XRD), respec- 
tively. The hardness of processed layers was investigated by a 


Fig. 11. “Frosted” area at the wall of the runner blade trailing edge due to outlet 
traveling bubble cavitation, shaded area C of Fig. 8. 


Fig. 12. Typical erosion at the wall of the runner hub due to interblades cavitation 
vortices, shaded area D of Fig. 8. 


Vicker hardness tester. The results indicated that there was 
different susceptible to plastic deformation caused by cavitation 
loading for different kind of laser processing (Fig. 18). 

Nicolet et al. [26] investigated is a 1 GW, four Francis turbines 
power plant. Francis turbines feature a cavitating vortex rope in 
the draft tube resulting from the swirling flow of the runner outlet. 
For hydraulic system modeling hyperbolic differential equations 
are 


dh a dQ 
at gA ox” (10) 
dh 1 0Q_, AlQ| Q=0 (11) 


dx © gA dt © 2gpA? 


where ‘a’ is the wave speed (m/s), ‘Q’ is the discharge (m?/s), ‘h’ is 
the piezometric head (m), ‘D’ is the pipe diameter (m), ‘A’ is the 
cross-section (m7), ‘g’ is the gravity (m/s) and ‘A’ is the friction loss 
coefficient. The unsteady pressure field related to the precession of 
the vortex rope induces plane wave propagating in the entire 
hydraulic system. The simulation results revealed the piping 
natural frequencies that were excited by the draft tube pressure 
source. In addition, the transfer function between the draft tube 
pressure source and the generator electromagnetic torque pointed 
out the risk of electrical power swing. However, the risk can be 
really evaluated only knowing the pressure excitations and the 
draft tube wave speed. If the first one can be obtained from scaled 
model testing, the latter has to be estimated, either experimentally 
from vortex rope photography, or in the future, by CFD. The 
methodology was a helpful tool for predicting the risks of 


0,4 


Fig. 13. Limits of cavitation development within the operating range of a Francis 
turbine. (1) Suction side leading edge cavitation limit; (2) pressure side leading edge 
cavitation limit; (3) interblade cavitation vortices limit; (4) discharge ring swirl 
cavitation limits. 
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Fig. 14. Typical eroded areas of a Kaplan runner. 


resonance at the early stages of pre-design or as a help for on site 
diagnostic purposes. 


4. Analytical investigations 


Muntean et al. [27] carried out a numerical investigation of the 
cavitating flow in a Francis turbine runner. First, the steady non- 
cavitating relative flow was computed in a runner interblade 
channel using the mixing interface approach. Second, the 
cavitation model was activated. Results for cavity shape and 
extent, as well as for the pressure distribution on the blade with 
and without cavity were presented and discussed. The equation of 
pressure coefficient (C,) was as given below: 


p- Pref 
C = 12 
p pEref ( ) 


where Pref is the reference pressure, p is the density and E is the 
specific energy. 

Qian et al. [28] carried out a three-dimensional unsteady 
multiphase flow simulation in the whole passage of Francis 
hydraulic turbine. The pressure pulsation was predicted and 
compared with experimental data at different positions in the draft 
tube, in front of runner, guide vanes and at the inlet of the spiral 
case. The relationship between pressure pulsation in the whole 
passage and air admission was analyzed. The computational 
results showed that air admission from spindle hole decreased the 
pressure difference in the horizontal section of draft tube, which in 
turn decreased the amplitude of low-frequency pressure pulsation 
in the draft tube; the rotor-stator interaction between the air inlet 


Fig. 15. Cavitation at the blade root. 
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Fig. 16. Relationship between inlet and outlet diameter of Francis turbine. 


and the runner increases the blade-frequency pressure pulsation in 
front of the runner. 

Bajic [29] worked out novel technique for vibro-acoustical 
diagnostics of turbine cavitation and its use demonstrated on a 
Francis turbine. The technique enabled identification of different 
cavitation mechanisms functioning in a turbine and delivered 
detailed turbine cavitation characteristics, for each of the 
mechanisms or for the total cavitation. The characteristics 
specified the contribution of every critical turbine part to the 
cavitation intensity. Typical diagnostic results were obtained as 
enabled optimization of turbine operation with respect to 
cavitation erosion, showed it how a turbine’s cavitation behavior 
can be improved; and formed the basis for setting up a high- 
sensitivity, reliable cavitation monitoring system. 

Yuliun et al. [30] carried out a study of cavitating flow in a 
Kaplan turbine having numerical simulation with a cavitation 
model and a mixture two-phase flow model, which are in 
commercial code Fluent 6.1. It was performed to unsteady 
cavitating turbulent flow within the entire flow passage of a 
Kaplan turbine in a hydraulic power plant. Based on the calculation 
results, region and degree of cavitation occurred in the flow 
passage of the turbine running in a given working conditions. 

Muntean et al. [31] carried out numerical investigation of the 
3D flow in the distributor (stay vanes and guide vanes) of the 
GAMM Francis turbine. The domain corresponded to the dis- 
tributor (stay vane and guide vane) interblade channel. There were 


60 
50 
z 40 
= Draft pipe 
ce k 
28 30 i 
g 
a —— w 
= xp w 
A > l . 
Main guide E 
10 + 


beaning case 
a 


pn — 
oo 5 10 15 20 25 30 85 
Output power 
(MW) 


=- : water-turbine runner with splitter blades = : existing runner 


Fig. 17. Results on vibration between turbine metal case and draft tube. 
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Fig. 18. Sample surface after grinding and polishing—(a) after 2 min, (b) after 4 min, (c) after 5 min and (d) after 7 min of cavitation erosion. 


three main issues addressed in the paper: first, using the numerical 
methodology presented above, the distributor flow for several guide 
vane angle values was computed as shown in Fig. 19. As a result, the 
guide vane torque versus guide vane opening angle was computed 
for the actual position of the guide vane axis; second, it was 
investigated the flow for the whole range of the guide vane positions, 
at four different locations of the guide vane axis. The optimization 
criterion considered was the minimization of the mechanical 
loading of the turbine regulating system. This means that the 
extreme value of the torque applied to the guide vane shaft has to be 
minimized. The designed parameter to be optimized is the location 
of the guide vane shaft axis along the guide vane chord line. 
Nilsson and Davidson [32] a parallel multiblock finite volume 
CFD code CALC-PMB (parallel multiblock) for computations 
of turbulent flow in complex domains was developed for an 
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Fig. 19. Three-dimensional computational domain for the GAMM Francis turbine 
distributor. 


investigation of the turbulent flow in Kaplan water turbines. This 
work was focused on tip clearance losses, which reduced 
the efficiency of a Kaplan water turbine by about 0.5%. The 
computational results from four different operating conditions, 
with different guide vane angles, were compared in the work. It 
was found that, the computations capture a vortical structure close 
to the leading edge tip clearance, where the tip clearance flow 
interacts with the shroud boundary layer and cavitational bubbles 
are formed. The tip blade loading increased when the specific 
speed decreased. Detailed measurements need to be performed, 
for use as boundary conditions and validation of the results. 

Gagnon and Deschenes [33] focused on unsteady rotor-stator 
interaction in a propeller axial turbine. The flow behavior 
was analyzed at different rotor and stator relative locations 
with numerical simulations using a commercial code and k-e 
turbulence model. The main goal was to study unsteady flow 
phenomena such as wake, separation, forces and pressure 
fluctuations in the propeller turbine. This investigation will help 
to design a series of flow measurements used in turn to improve 
future CFD simulations with realistic velocity profiles as boundary 
conditions. Wakes and unsteady rotor-stator interactions for a 
different operating regime of a propeller turbine was studied. The 
wake behind guide vanes was dissipating very fast. It was not 
predicted whether it had an effect on the rotor stator unsteady 
interactions since higher harmonics of the spectrum were not 
well defined. The main interaction was therefore attributed to the 
pressure field fluctuations. Finally, the upcoming experimental 
phase to increase our knowledge and understanding of the 
propeller turbine flow may be looked forward. 

Susan-Resiga et al. [34] carried out a numerical analysis of the 
swirling flow downstream through a Francis turbine runner. It was 
found that the flow stability characteristics were changed when 
decreasing the discharge. In Fig. 20 it was shown that the swirling 
flow in the survey section downstream the runner, in the draft tube 
cone, reached a critical state in the neighborhood of the best 
efficiency operating point. For larger discharge, the swirling flow 
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Fig. 20. Jet control technique for swirling flow in the discharge cone of Francis 
turbines. 


was supercritical, and thus it was not able to sustain ax 
symmetrical perturbations. However, at partial discharge the flow 
became subcritical and it was able to sustain ax symmetric 
perturbations. Further investigations revealed that the axial 
velocity and specific energy deficit in the central region were 
responsible for the helical vortex breakdown, also known as 
“precessing vortex rope”. 

Aschenbrenner et al. [35] used the mixture model and full 
cavitation model to compute the unsteady turbulent flow and 
compared with experiment. From the vaporization and condensa- 
tion mechanics between vapor and liquid and the assumption of 
unchanged pressure in cavity, the mass change rate of cavity was 
proportional to liquid pressure and vapor pressure difference. 

Nennemann et al. [36] studied blade surface and discharge ring 
cavitation with an emphasis on the latter using CFD. Localized high 
flow gradients (cavitation) as well as more global unsteady effects 
(rotor stator interaction) contributed to the phenomena requiring 
advanced and resource intensive CFD approaches, notably locally 
refined meshes coupled with unsteady computations including 
runner rotation. CFD results enabled to explain the presence of 
discharge ring cavitation at distinct circumferential locations 
corresponding to the number of guide vanes used unsteady 


Fig. 21. Tip gap and tip gap vortex cavitation on a model Kaplan runner. 
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Fig. 22. Cavitation details in the runner domain in high cavitating operation. 


pressure measurements and high-speed camera images to 
qualitatively support the validity of our CFD model. Ultimately 
it was found that the CFD results enabled to understand the physics 
underlying the problem of discharge ring cavitation. Extensive 
parameters were studied to investigate all the significant influence 
factors as shown in Fig. 21. 

Balint [37] carried out a numerical investigation by computing 
the 3D turbulent single phase flow in the Kaplan turbine runner, 
the secondary phase for the water (vapor) was activated starting 
far from the cavitation occurrence conditions. The output pressure 
in the runner domain was decreased to initiate the cavitating 
conditions on the runner blades. The unsteady solver with reduced 
time step size was employed to capture the self-induced cavitation 
in the runner domain and 3D effects together with wall-friction 
modification onto the runner blades. The shock of the unsteady 
cavitating flow leads to high distortion of the flow ingested 
downstream by the draft tube of the turbine. It was concluded that 
unsteady effects of the flow have been made mainly by the 
unsteady detachments of the cavitation at the blade suction side 
close to the trailing edge. An original method of computing the 
cavitation inception and of comparing different operating regimes 
in similar cavitating conditions was presented by plotting the 
cavitation volume growth as shown in Fig. 22. 


5. Case studies 
Bauer et al. [38] explored the techniques for the purpose of 


visualizing isolated flow structures in time-dependent data. 
Primary industrial application was the visualization of the vortex 


Fig. 23. Cavitation bubbles near Kaplan runner blades. 
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Fig. 24. Vortex rope in Francis draft tube. 


rope, a rotating helical structure which builds up in the draft tube 
of a water turbine. The vortex rope can be characterized by high 
values of normalized helicity, which is a scalar field derived from 
the given CFD velocity data. In two related applications, the goal 
was to visualize the cavitation regions near the runner blades of a 
Kaplan turbine and a water pump as shown in Fig. 23. Again, the 
flow structure of interest can be defined by a scalar field, namely by 
low pressure values. A particle seeding scheme based on quasi- 
random numbers was proposed, which minimized visual artifacts 
such as clusters or patterns. By constraining the visualization to a 
region of interest, problems were reduced and storage efficiency 
was gained. They applied method to various data sets from our 
industry partners, visualizing the vortex rope in the draft tube of a 
Francis turbine shown in Fig. 24, and cavitation on the suction side 
of various turbine and pump runner blades. The concept of 
selectively visualizing relevant flow structures was proven to give 
additional insight into their complicated dynamic behavior. 


6. Conclusions 


Cavitation is a phenomenon of formation of vapor bubbles in 
low pressure regions and collapse in high pressure regions, high 
pressure is produced and metallic surfaces are subjected to high 
local stresses. It is difficult to avoid cavitation in hydro turbines 
which cannot be avoided completely, but can be reduced to an 
economically acceptable level. Many investigators have studied 
the process of cavitation in hydro turbines through experimental 
and analytical studies. Some of the investigators have reported that 
in spite of design changes in the turbine components and providing 
different materials and coatings to the turbine blades, however the 
improvement in most cases is not quite significant. It is therefore; 
further experimental and theoretical studies are required for 
studying the impact of cavitation in hydro turbine to determine 
impact of cavitation at different values of parameters which was 
relate to the cavitation in hydro turbine. CFD based analysis of 
cavitation in reaction turbines could be cost effective solution for 
an extensive analysis. 
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